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Abstract 

We study the motion of vortices in the conserved and non-conserved phase- 
ordering models. We give an analytical method for computing the speed and 
position distribution functions for pairs of annihilating point vortices based on 
heuristic scaling arguments. In the non-conserved case this method produces 
a speed distribution function consistent with previous analytic results. As two 
special examples, we simulate the conserved and non-conserved 0(2) model 
in two dimensional space numerically. The numerical results for the non- 
conserved case are consistent with the theoretical predictions. The speed 
distribution of the vortices in the conserved case is measured for the first 
time. Our theory produces a distribution function with the correct large speed 
tail but does not accurately describe the numerical data at small speeds. The 
position distribution functions for both models are measured for the first time 
and we find good agreement with our analytic results. We are also able to 
extend this method to models with a scalar order parameter. 
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I. INTRODUCTION 



The phase ordering dynamics of certain physical systems after a rapid temperature 
quench below its critical temperature is dominated by the annihilation of topological de- 
fects of opposite charge [1]. In particular the n- vector model with non-conserved order 
parameter (NCOP) Langevin dynamics, where the defects are vortices, has been studied in 
some detail [2-4]. 

Mazenko [4,5] carried out an investigation of the distribution of defect velocities for non- 
conserved phase ordering systems. By using an approximate "Gaussian closure" scheme, 
he was able to compute the velocity distribution for vortices in the non-conserved n-vector 
Langevin model for the case of point defects where n = d dimensions. We [7] carried out 
numerical simulations for the n = d = 2 non-conserved case and measured the vortices 
speed distribution. The results are consistent with Mazenko's theoretical predictions. In 
particular the power-law tail of the distribution at large speeds which is robust is correctly 
predicted. The problem of the relative velocity as a function of separations for annihilating 
pairs was treated in Ref. [5], and the velocity distribution for strings for the non-conserved 
order parameter case was treated in Ref. [6]. 

Bray [8] developed a heuristic scaling treatment of the large speed tails based on the 
disappearance of small defects (annihilating pairs or contracting compact domains). This 
method treated only the power-law exponent of the distribution's large speed tail. For the 
non-conserved n-vector case this simple argument gives a result consistent with Mazenko's 
theory. However this method is also able to produce the large speed tail exponents for the 
conserved n-vector models, and conserved and non-conserved scalar (n=l) models. 

We show here that Bray's arguments can be extended to give results beyond the tail 
exponents. In models where point defects dominate the dynamics, one can compute the 
defect speed distribution functions based on Bray's scaling assumption. The same idea can 
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be easily generalized to the scalar order parameter case. 

In a very recent work Mazenko [9] has suggested how the previous work in Ref. [4] can 
be extended to anisotropic systems and the conserved order parameter (COP) case. He 
finds that the average speed goes as t^ 1 for the COP case with a scaling function of the 
same formal form as for the NCOP case given by Eq. (40) below. These results are not 
in agreement with the analytical or numerical work presented below in this paper. The 
Gaussian closure method developed in Ref. [9] does not appear adequate for treating the 
COP case. 

In the next section we will generalize Bray's argument for the point defect case. We find 
simple analytic expressions for the speed and separation distribution functions. We recover 
the large speed tail exponents obtained previously by Bray. For the non-conserved n-vector 
model we obtain precisely the same results as found in Ref. [4]. Then in Sec. Ill we present 
the numerical simulation results for non-conserved n = d = 2 Langevin model. Next in Sec. 

IV, we present the simulations results for the conserved n = d = 2 Langevin model. In Sec. 

V, we point out that the method developed in this paper can be used for those cases where 
one has a scalar order parameter. 



II. THEORETICAL DEVELOPMENT 

Let us suppose that we have N pairs of oppositely charged vortices which are on their 
way to annihilation. We suppose that pair % is separated by distance rj(t) with relative speed 
Vi(t) = \fi(t)\. Consider the associated phase-space distribution function: 

N 

f(r,v,t) = (£ l 8(r-r i (t))6(v-v i (t))) . (1) 
i=i 

This quantity satisfies the equation of motion 

r) r) N r) N 

m f(r, v, t) = --(£ mS(r - ri (t))Hv - «*(*))> -^<E Vi(t)Hr ~ n(t))5(v - Vi (t))) . (2) 
i=i i=i 



Our key kinematical assumption is that the relative velocity is a known function of the 
separation: 

hit) = - Vi (t) = -u(n(t)) (3) 



ii(t) = u'inityfiit) = -u(ri(t))u'(ri(t)) . (4) 

We check these assumptions as we proceed. Eq.(2) then takes the form 


—f(r,v,t) = —(u(r)f(r,v,t)) + — (u(r)u'(r)f(r,v,t)) . (5) 
where we have the normalization 



oo poo 



dr / dvf(r,v,t) = N(t) . (6) 
o Jo 

Eq.(5) is one of our primary results. 

Our assumptions are consistent with being in a regime where the annihilating pairs are 
independent and we can write 

f(r,v,t) = N(t)P(r,v,t) (7) 

where P(r,v,t) has the interpretation as the probability that at time t we have a pair 
separated by a distance r with relative speed v. Inserting Eq.(7) into Eq.(5) we find that 
P(r, v, t) satisfies 

—P(r, v,t) = — (u(r)P(r, v, t)) + — (u(r)u'(r)P(r, v, t)) + 7 P(r, v, t) (8) 

where 

i=~m) ■ ( 9 ) 

We will see that 7 (and N(t)) are determined self-consistently by using scaling ideas. 
We are interested in the reduced probability distributions 



poo 

P r (r,t) = / dvP{r,v,t) , (10) 
Jo 

and 

POO 

P v (v,t) = / drP(r,v,t) (11) 
Jo 

with the overall normalization 

POO POO 

dr dvP(r,v,t) = 1 . (12) 
Jo Jo 

Our goal is to solve Eq.(8). The first step is to show that 

P(r,v,t) = P r (r,t)8(v -u(r)) . (13) 

Inserting Eq.(13) into Eq.(8) we have 

d d 
S( v - u(r))—P r (r, t) = ~f5(v - u(r))P r (r, t) + — (u(r)5(v - u(r))P r (r, t)) 

d 

+— (u(r)u'(r)5(v — u(r))P r (r, t)) 

= 6(v-u(r)) (jP r (r,t) + ^ (u(r)P r (r, t))) 

+u(r)P r (r, t) (J^S{v - u(r)) + -^u'(r)5(v - u(r)^j . (14) 
Using the following identity 

0-5(v-u(r)) = -?- [ -e iX ^)) 
dr v v " dr J 2vr 

J 27T 

= -u'(r)^-6(v-u(r)) (15) 
ov 



we find that Eq.(13) holds with P r (r,t) determined by 

^P r (r, t) = 7 P r (r, t) + ^ (u(r)P r (r, t)) . (16) 
Imposing the normalization, 

poo 

/ drP r (r,t) = l , (17) 
Jo 

we find on integrating Eq. (16) over r that 

7 = ]hnu(r)P r (r,t) . (18) 

r— >0 

Thus 7 and iV(£) are determined self-consistently in terms of the solution to Eq. (16). 

So far this has been for general w(r), let us restrict our subsequent work to the class of 
models where the relative velocity is a power law in the separation distance: 

u = Ar-- b (19) 

where A and b are positive. Next we assume that we can find a scaling solution [11] to Eq. 
(16) of the form 

Pr(r,t) = j^F(r/L(t)) (20) 

where the growth law L{t) is to be determined. Inserting this ansatz into Eq.(16) we obtain: 

—L b L (xF' + F) = Ax- b F' - Abx~ b ~ l F + L 6+1 7 F (21) 
where x = r/L(t). To achieve a scaling solution we require 

L b L = AC (22) 

and 

L b+1 -f{t) = AD (23) 
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where C and D are time independent positive constants, the factors of A are included for 
convenience. Eq. (21) then takes the form 

-C [xF' + F}=DF + x~ b F' - bx-^F . (24) 

This has a solution 

F <*> = (irl^F (25) 

where B is an over all positive constant and the exponent in the denominator is given by 

a = l + a/z (26) 

where z = 1 + b and a = D/C. If we enforce the normalization Eq. (17), we find that 
D = B. This reduces the spatial probability distribution to a function of two unknown 
parameters B and C assuming that b is known. 

We are at the stage where we can determine the number of annihilating vortex pairs as 
a function of time. From Eqs.(18), (9) and (25) we have that 

AB AD 

7(0 = JT+b = JT+b (27) 
and we again have that D = B. However from Eq.(22) we have 

J_± L i+b = A c (28) 

1 + bdt K J 

and for long times 

L 1+b = ACzt . (29) 

Putting this result back into Eq.(27) we find 

, , AB a . . 

^ = Acx = 7t ■ (30) 

We have then, from Eq. (9), that the number of pairs of vortices as a function of time is 
given by 



N(t)/N(t ) = (t /t) a/z . (31) 

However from simple scaling ideas we have rather generally that for a set of point defects in 
d dimensions 

N(t) « L- d . (32) 
Comparing this with Eq.(31) we identify a = d. This gives our final form for F(x) 

F ^ = (i + c%)^/' ■ (33) 

We check the validity of this result in sections III and IV. 
The speed probability distribution is given by 



oo 



P v (v,t) = / dr P(r,v,t) 



o 



POO 

/ dr S(v- Ar- b )P r (r,t) 
Jo 



B 1 / C , 

1 + n^T^m ( 34 ) 



AbL 1+b v 2+1 / b \ L 1 + b v( 1+b y b 
where v = v/A. We can define the characteristic speed v via 



,£v (l+6)/6 _ C 

- L l+^(l+6)/& ^5) 



or 

C(l+6)/6 



and 



«(*) = — « *~ (1 ~ 1/J ° ( 36 ) 



^M) = ( X + = Jp.(t,/«) (37) 
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where V — v/v and the distribution function has a scaling form. Clearly the large speed 
tail goes as V~ p where p = 2 + 1/6 in agreement with Bray's result. After rearrangement 
we find 

PvM = l Pv{v/f}) = lL (i + vw>y (38) 

where 

• = i-l = f-l • (39) 
We numerically test this result for various models below. 

III. NON-CONSERVED N- VECTOR MODEL 

We now want to test our theoretical results for P r {r, t) and P v {v, t) for the non-conserved 
time-dependent Ginzburg-Landau (TDGL) 0(n) model where b = z — 1 = 1. If we work in 
terms of dimensionless variables v — v/v and r — r/f where v(t) and f(t) are the average 
speed and separation as function of time, then Eq. (38) gives the vortex speed distribution 
function 

P.(5) = »3" /2 (1+ ^„ +2)/2 • («) 
with (5 = 7r[r(iin)/r(|)] 2 . This is exactly the familiar result found in Ref. [4] for n = d. 
The average speed is v oc t^ 1 / 2 and z = b + 1 = 2. 
As a special case, when n = 2, we have 

(1 + pv 2 ) 

with (3 = (vr/2) 2 = 2.4674. Both p v (v) and v(t) have been verified in Ref. [7]. The energy 
and defect number are proportional to (t/lnt) -1 , where there is a logarithmic correction. 
But we did not see such a correction for the average speed v(t) oc t" 1 ^ 2 . 
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Let us turn next to the distance distribution function. From Eq. (33) we have, for 
n = d = 2, 

ICr 

F ^ = (i + cP? ' (42) 

where C = 2.4674. 

We check this numerically using the same data as in Ref. [7]. The model is described by 
a Langevin equation defined in a two-dimensional space 

— * 

| = e^+vV-WV, (43) 

where e is set to be 0.1, and the quench is to zero temperature, so we need not include 
noise. We worked on 1024 x 1024 system with lattice spacing Ar = 7r/4. Periodic boundary 
conditions are used. Starting from a completely disordered state, we used the Euler method 
to drive the system to evolve in time with time step At = 0.02. 

The position of a vortex is given by the center of its core region, which is the set of points 
(xi,yi) that satisfy \ip(xi,yi)\ < (H)/4. By fitting \ip(xi,yi)\, where {x h yi) are the points 
belonging to a vortex's core region, to the function M(x, y) = A + B[(x — x ) 2 + (y — yo) 2 } we 
can find the center (x , y ). The positions of each vortex at different times are recorded, and 
the speed is calculated using v = Ad/ At. Here Ad is the distance that the vortex travels 
in time Ar = 5. 

To measure f{t) and F(r) we must first accumulate the following data. In a given run 
we keep track of the trajectories of all the vortex centers. We label each pair of oppositely 
charged vortices which annihilate. Then move backwards in time to determine for each such 
pair the separation as a function of time Ti(t). Then r(t) is the average separation between 
annihilating pair of vortices at time t. The average distance f(t) is shown in Fig. 1. From 
the discussion in Sec. II we expect 

f(t)=L(t) = a[t+(f /aY] 1/z , (44) 
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where r = r(t = 0) and a is a constant. The average distance between annihilating pairs 
increases with time and a fit to the data gives a = 2.71, f 50 and z ~ 2.22. 
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FIGURES 



100 
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FIG. 1. The average distance f between annihilating pairs versus time after quench. The data 
is averaged over 68 runs. The fit to the data is given by Eq. (44). 

To measure the probability distribution function F(r), we distribute the various pairs 
into bins of width A = 0.01 centered about the scaled separation ri(t)/f(t). We then plot 
the number of pairs in each bin versus r(t)/f(t) and properly normalize to obtain the scaling 
result shown in Fig. 2. In the following, when we measure the other distribution functions 
with scaling properties we employ the same method. The curve representing F(x) given by 
Eq. (42) is also shown in Fig. 2. There is no free parameter in the fit other than b and 
a. The fit is fairly good. At large distances we can see the function approximately obeys a 
power law. The exponent is about 6, which is different from the value 3 indicated by Eq. 
(42). We do not know why it is so. 
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FIG. 2. Separation probability distribution F(r) versus the scaled separation for the NCOP 
case with n = d = 2. The data is averaged over 68 runs with a bin size of 0.01. The solid line is 
Eq. (42) with 6=1 and a = 2. In the insert we show the same data on a logarithmic scale. At 
large r the distribution is approximately a power law with an exponent about 6. The dashed line 
in the insert is proportional to r~ 6 . 

We next measure u(r), the average speed for annihilating defects separating by a distance 
r. We track the motion of each annihilating pair, and determine for each pair the speed 
Ui = (ri(t + At) — r,i(t))/ At as a function of r^. Then we average Ui(r) over all the pairs 
that have a fixed r\ = r. The result is shown in Fig. 3. For small enough separations we 
have u(r) oc r~ b where b ~ 1 as expected. 
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FIG. 3. Test of Eq. (19). See text for a discussion. The data is averaged over 68 runs. There 
is a scaling regime at small distances r with an exponent near to —1. At large distances, while the 
statistics are not as good, there is still an approximately power-law dependence on r. 



IV. CONSERVED N- VECTOR MODEL 

Let us turn to the case of a conserved order parameter where for a TDGL model we 
expect b = z — 1 = 3. In this case the vortex speed distribution, Eq. (38), is given by 

,~ nf3 n ' A #»-3)/3 

MV) = — (l + /W 3)<^>/< (45) 

with f3 = (nr(|)r( 3 -^)/r(l + |)) 4/3 and the average speed is v ~ t~ 3/4 . As a special case, 
consider n = 2, where 

... g-' /3 , ... 

* W = ~ (1 + ,«W>)"" <46) 
with /9 = 2.27773.... Notice, unlike the NCOP case, p v (v) blows up for small v. This appears 

to be an unphysical feature. 

The distribution function for the distance between annihilating pairs, Eq. (33) with 

b = 3 and n = d = 2, gives 

F « = (TT§^ • < 47 > 

where C = 11.817.... 

We simulated the conserved 0(2) model in two dimensions to test these predictions. The 
model is described by a Langevin equation denned on a two-dimensional space 



at L 

= V 2 ^M , (48) 
5^ 



where the effective Hamiltonian is given by He[ 



d 2 r. 



All the quantities are dimensionless. We work on 256 x 256 system with lattice spacing 
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Ar = 7r/4 and again periodic boundary conditions are used. We employ the method invented 
by Vollmayr-Lee and Rutenberg [10] to numerically integrate Eq. (48). This method is stable 
for any value of integration time step At. As the time t increases, the evolution of the system 
becomes progressively slower. With the new time step technique we can increase the time 
step to accelerate the evolution. We let At = 0.01 t a36 after t > 120. 

In addition to the vortex statistics discussed in the NCOP case, we also measure the 
average energy E = (He) above the ground state energy E = —S/A with S being the area 
of the system. The energy and number of vortices are shown in Fig. 4. The power-law 
exponent for the defect number is —0.51, which is consistent with d/z = 1/2. The decay 
power law for the energy is t~ 0A1 , which is slower than that for the defect number. This 
may be due to the relaxation of spin waves. 
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FIG. 4. The energy E(t) — Eo and the vortex number N(t) are plotted versus time after quench. 
The data is averaged over 61 runs. The dashed lines, which are guides to the eye, are proportional 
to t~ 0A1 for the energy plot and i -0 - 51 for the vortex number plot respectively. 

We use the same method as in the non-conserved case to find the center for each vortex. 
The speed of each vortex is computed by using v = Ad/ At with r = 10. We measure the 
speed for each vortex at the same time t and average over different vortices. The average 
speed of the defects is shown in Fig. 5. The prediction for the exponent is —(1 — 1/z) = 
—0.75, while the measurement finds —0.77. 

o.i 

o.oi 

i > 

0.001 

aooo roo 1000 10000 

t 

FIG. 5. The averaged speed of vortices at time t for the COP case. The data is averaged over 
61 runs. 

Next we determine the speed distribution as a function of time. Again we plot the scaled 
data from different times to test the scaling property of this distribution function. The 
resultant p v (v) is shown in Fig. 6. Clearly scaling works and the large speed tail exponent is 
2.24. This is close to the prediction 2 + 1/b = 2.33. However the theory fails at small scaled 
speeds where the simulations go to zero while the theory blows up. Clearly the exponent s 
in Eq. (38) is poorly determined in the theory for b = 2 and d = 2. If we allow b and a float 
then we obtain an excellent fit shown in Fig. 6. 
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FIG. 6. The probability distribution for vortex speed for the COP case. The data are from 
61 different runs. The bin size is 0.01. The solid line is Eq. (46). The dashed line is the best fit 
to p v (v) by changing b and a. In the insert we show the same data on a logarithmic scale. The 
dashed line in the insert is used to guide the eye and proportional to v~ 2 - 24 . 

The average separation of annihilating pairs of vortices f{t) for the COP case is shown 
in Fig. 7. Again fitting this to the form given by Eq. (44) we find a = 2.285, f = 17 and 
z = 4.0. 
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FIG. 7. The average distance r between annihilating pair of vortices versus time after quench 
for the COP case. The data is from 61 runs. The fit to Eq. (44) is given in the text. 
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As in the NCOP case we can measure the separation distribution function F(r). This is 
shown in Fig. 8. With no free parameters, b and a being fixed, the fit is pretty good. At 
large distances, the statistics are poor. But we can see the function approximately obeys a 
power law. The exponent is about 4, which is different from the value 3 indicated by Eq. 



FIG. 8. The probability distribution for the separation between two annihilating vortices for 
the COP case. The data is averaged over 61 runs. Bin size is 0.01. The solid line is Eq. (47) with 
6 = 3 and a = 2. In the insert we show the same data with a logarithmic scale. The behavior 
of this function at large distances is approximately a power law with an exponent about 4. The 
dashed line is proportional to f -4 . 

We also measure u(r), as in the NCOP case. Our results are shown in Fig. 9. The 
assumption u ~ r~ h is well satisfied with b m 3 at small enough distances. 



(47). 




r / r 



18 




FIG. 9. The average speed as a function of pair separation r for the COP case. There is a 
scaling regime at small r with the exponent near to b = 3. The data is collected over 61 runs. 

V. SCALAR MODELS 



We follow the argument given by Bray [8] to extend our discussion to include models 
with a scalar order parameter (n = 1). First we calculate the probability function P r (r,t) 
for the domains with radius r. This calculation is the same with that in Sec. II. We obtain 
Eq. (33) for F(x). Next, we compute the area- weighted probability for interfacial radius of 
curvature r by multiplying F(x) with x d_1 and then normalize the resulting quantity. The 
resulting probability function F s (x) is 



F s (x) 



dr(d/z)c( 1+d+z y z 



X 



d+z-2 



r(i/2)r((d + z - i)i ' z ) (i + Cx z y+ d / z ' 

Following Bray we use v = Ar b to get the distribution function for the interface speed 



(49) 



Pv(v) = 



dr(d/z)pV z 



v 



-1 + 1/(2-1) 



(50) 



(z - i)r(i/*)r((d + z - 1)/ 'z) (i + /3£fi+i/(*-i))i-M/* • 

For the non-conserved case, this result is the same with the one obtained by Bray using the 
Gaussian calculation. The large speed tail exponent is p = 2 + d/(z — 1), which is valid for 
both conserved and non-conserved models. 
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VI. CONCLUSION 



We show how a simple generalization of Bray's scaling argument can lead to quantitative 
results for certain distribution functions. In particular we find that the distribution function 
for the distance between annihilating pairs of vortices is well described by the scaling theory 
for both NCOP and COP dynamics for n — d — 2. We are also able to compute the 
speed distribution function using these ideas. For non-conserved models, we reproduce the 
accurate result obtained previously. For conserved models, the speed distribution function 
only gives us the correct tail exponent. 

Our method can also be extended to scalar cases, and generate a full expression for the 
interfacial speed distribution. The power-law tail exponent p = 2 + d/(z — 1) is obtained. 
The result is the same as the result obtained by Bray [8] . 

The simple scaling method presented here leads to a reasonable description of the statis- 
tics of defect dynamics. Clearly it does a better job for the NCOP case since the speed 
distribution function for COP case does not show the proper small speed behavior. Simi- 
larly the more microscopic method of Ref. [9] leads to an adequate treatment of the small 
speed regime for the COP case but does not give the correct large speed tail. 

The approach developed here is highly heuristic. Can it be systematized? Clearly to 
improve this approach one would need to include the interactions between different pairs. 
It is not clear how one does this. 
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